Non-Abelian quantized Hall states of electrons at filling factors 12/5 and 13/5 

in the first excited Landau level 
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We present results of extensive numerical calculations on the ground state of electrons in the first 
excited (n — 1) Landau level with Coulomb interactions, and including non-zero thickness effects, for 
filling factors 12/5 and 13/5 in the torus geometry. In a region that includes these experimentally- 
relevant values, we find that the energy spectrum and the overlaps with the trial states support the 
previous hypothesis that the system is in the non-Abelian k = 3 liquid phase we introduced in a 
previous paper. 
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Many distinct quantum Hall liquid phases have been 
observed in two-dimensional electron systems. Leaving 
aside those that occur at integer filling factor (or quan- 
tized Hall conductance), the non- integer fractions occur 
at low filling factors in high mobility samples (at higher 
filling factors, they are supplanted by states in which 
translational or rotational symmetry is violated, or by 
"re-entrant" integer quantized Hall phases). In the low- 
est (n — 0) Landau level (LL), that is at filling factors less 
than 2, the incompressible liquids are phases of matter 
that are correctly characterized by the Laughlin states 
Q , and their extensions via the hierarchy Q or compos- 
ite fermion 0] approaches; these two approaches describe 
the same phases jj|. In the first excited (n = 1) Landau 
level, the physics appears to be different, as one gets 
closer to the broken symmetry phases. In general, frac- 
tions in the n = 1 LL can be compared with those in the 
n = level by subtracting 2 from the filling factor; this 
corresponds to simply filling the lowest level with both 
spins, and treating the next level like the lowest. A quan- 
tized Hall plateau occurs at fillings 5/2 and 7/2, which 
is now widely believed [5, 6] to correspond to an incom- 
pressible liquid that may be viewed as a p-wave paired 
state of spin-polarized composite fermions in a half-filled 
LL 0. It is unlike the lowest Landau level (LLL) case, 
in which at filling factor v = 1/2 or 3/2 the composite 
fermion liquid appears so far to be gapless, and exhibits 
Fermi-liquid-like properties. 

A few other fractions are observed between 2 and 4, 
and one may wonder whether these are the same phases 
as occur in the LLL. In an earlier paper |a] (to be re- 
ferred to as RR), a sequence of incompressible fractional 
quantum Hall liquids was constructed. These have filling 
factors v — k/(Mk + 2), where fc = 1, 2, 3, and 
M = 0, 1, 2, are integers, and for fermions M must be 
odd (the values [(M - l)k + 2}/(Mk + 2) may also be 
obtained, by using particle- hole symmetry). These were 
constructed within the LLL, but may be applied to higher 
filling factors by adding the filling of the lower levels. The 
k = 1 liquids are the familiar Laughlin states Q, while 
fc = 2 is the Moore- Read (MR) paired state 0. For the 
next case, k = 3, some evidence that it occurs in the 
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was presented, but was perhaps not entirely convincing. 
Meanwhile, experiments have observed a fractional quan- 
tum Hall state at 12/5, which exhibits a remarkably small 
energy gap for charged excitations 0- There is a great 
deal of interest in phases like the RR states, as for fc > 1 
they exhibit excitations that obey non-Abelian statistics 
[7j, and for fc ^ 1, 2 or 4 they would support universal 
quantum computation 

In this paper we return to the issue of the nature of 
the liquid state at filling fractions 12/5 and 13/5. We 
report numerical calculations for moderate numbers of 
electrons in a single LL, with interactions that model 
the Coulomb interaction between electrons in the n = 1 
LL, plus the effect of non-zero thickness of the electron 
wavefunction in the direction perpendicular to the two- 
dimensional layer. We also explore the phase diagram as 
the short-range component of the interaction is varied. 
For experimentally relevant values of thickness, we find 
that the liquid phase in the vicinity of the n = 1 Coulomb 
interaction appears to be the RR fc = 3 phase. The evi- 
dence for this comes from the spectrum on the torus (i.e. 
periodic boundary conditions on a parallelogram), which 
exhibits a doublet of ground states that is characteristic 
of the fc = 3 phase, separated from higher excited states 
by a significant gap; this doublet does not occur for the 
hierarchy/compositc-fcrmion (H-CF) phase. Quantita- 
tively, the gap in the spectrum is small, suggesting that 
the gap for charged excitations will also be much smaller 
than for the H-CF 2/5 state in the LLL, in general agree- 
ment with experiment. The states in the ground state 
doublet have large overlaps with the trial states of RR, 
adapted to the torus. When the interaction is modified, 
phase transitions to a "stripe" phase, and to the usual 
H-CF phase are observed. Several of these results, except 
possibly for this last transition, are similar to results of 
a recent study of bosons in the LLL at filling 3/2, 
which is the k = 3, M = case of RR, and also to the 
1/2 (or 5/2, etc) case for electrons Q. 

The methods of calculation are rather standard, so we 
will be brief. For spin-polarized electrons in two dimen- 
sions confined to any one LL, the interaction Hamiltonian 
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i?i nt in the infinite plane can be represented by the pseu- 
dopotentials V m , rn = 1, 3, .... V m is the interaction 
energy for a single pair of electrons of relative angular 
momentum to (to odd) |l2j . For the case of zero thick- 
ness, the pseudopotentials for the Coulomb interaction in 
the n = 1 LL were plotted in Ref. 0. We also include 
non-zero thickness through the standard Fang-Howard 
method, in which the wavefunction contains the depen- 
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(but vanishes for £3 < 0) on the 



perpendicular coordinate X3 . This tends to reduce some- 
what the low to pseudopotentials relative to the others. 
We can then represent any one LL using the states in the 
LLL. This description of the interaction can be extended 
to the sphere (also using rotation symmetry), and to the 
torus; in these cases there are flux quanta piercing the 
system. Starting from the pseudopotentials for the non- 
zero thickness n = 1 LL, we also consider changing the 
first two pseudopotentials Vi, V3 by small amounts SVi, 
6V3 in order to explore the phase diagram in the vicinity 
of this interaction. This mitigates to some extent our ig- 
norance of the precise interaction Hamiltonian. For com- 
parison, for 15 particles on the torus and w = 2£b, the 
unperturbed values are V\ — 0.3858, V3 = 0.3333. We 
also note that particle-hole symmetry holds as long as 
inter-LL interactions are neglected, as here, so that our 
results for v = 13/5 also apply to 12/5, with no modifi- 
cations at all in the case of the torus geometry. Likewise, 
our results should also be relevant for filling factors 3 + | 
and 3 + |, which also lie in the n = 1 LL. Accordingly, 
we refer only to v = 3/5 from here on. 

We also refer to a (positive) p-body interaction which 
penalizes the closest approach of p fermions that is al- 
lowed by Fermi statistics; it can be written in terms of 
derivatives of (5-functions on the sphere or torus. For 
fermions (electrons), the parafermion trial states found 
in ref. Q are unique, exact zero-energy eigenstates of 
such interactions for p — k + 1 when N is divisible by 
k and = (k + 2)N/k — 3 (on the sphere), so that 
v = lim^oo N/N v = k/(k + 2) (thus M = 1). There are 
corresponding states on the torus, for Nj, = (k + 2)N/k. 
For the k + 1-body interaction, the trial states represent 
incompressible liquid phases, in which the excitations en- 
joy non-Abelian statistics for k > 1. The k + 1-body 
Hamiltonians allow us to numerically generate the trial 
states, for comparison with the exact ground states of the 
two-body pseudopotential interaction at the same N, 
in the same geometry. 

On the torus at v = N/N^ = 3/5, translational sym- 
metry implies that all energy eigenstates possess a triv- 
ial center-of-mass degeneracy of 5, which is exact for any 
size system, and also that there is a quantum number 
called K |b|, which is a vector lying in a certain Bril- 
louin zone. For K = 0, energy eigenstates can also be 
labeled by their eigenvalues under a rotation. In the 
M = 1 RR phases, the ground states have a net degen- 
eracy \{k + l)(fc + 2) in the thermodynamic limit, which 
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FIG. 1: Low-lying spectrum for 18 electrons on the torus 
on the torus vs. K = |K|, for hexagonal unit cell. Energies 
are in units of e 2 /£b- The interaction Hamiltonian is the 
n = 1 Coulomb interaction, including non-zero thickness with 
parameter w = 2£b, and SVi = —0.0035. 



is connected with the non-Abelian statistics [jj |£j • For 
k = 3, this 10- fold degeneracy is made up of the trivial 
factor 5 (which is always discarded in numerical stud- 
ies), together with a further 2-fold degeneracy, which in 
general becomes exact only in the thermodynamic limit; 
all these ground states have K = 0. (For the 4-body in- 
teraction, the 10-fold degeneracy is exact for any size, as 
the ground states have exactly zero energy.) By contrast, 
the H-CF ground state for fermions at v = 3/5 possesses 
only the 5-fold degeneracy. Then for an incompressible 
fluid on the torus, the spectrum of a sufficiently large 
system in one of these two phases should exhibit a nearly 
degenerate pair of ground-states or a single ground state, 
respectively, at K = 0, separated by a clear gap from a 
region of many states at higher energy eigenvalues. 

In Fig. ^ we show the spectrum for 18 electrons on 
the torus at 30 flux quanta, with thickness parameter 
w = 2b = 2, and a small change to the pseudopoten- 
tials, SVi = —0.0035, 5Vz = 0. In this and all subsequent 
figures, the geometry of the torus is the hexagonal unit 
cell. Energies are in units of e 2 /£b for the Coulomb in- 
teraction e 2 /r. We set the magnetic length £b and h to 
one throughout. For this ground-state doublet at 

K = is apparent, separated from the excited states by 
a gap which is about ten times larger than the splitting of 
the doublet, or than the typical spacings above the gap. 
This is the expected behavior of the RR incompressible 
fluid phase. However, we will discuss below the variation 
of this splitting with the interaction parameters. 

In Fig. we show the sum of the squared-overlaps of 
the lowest two K = states with the low-lying doublet in 
Fig. n with the two-dimensional subspace of zero-energy 
states of H4 (the trial states) for N — 15 particles, as 
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FIG. 2: The sum of the squared-overlap of each member of 
the ground-state doublet for n = 1 Coulomb interaction (for 
w = 2 and w) with the trial two-dimensional subspace, as a 
function of 8Vi on the torus, for 18 particles. [Also included 
is the same for a quantum well wavefunction in xz, labeled 
QW, with w = 2.75.] Also, the squared-overlap of one of 
the ground state doublet for w = 2 with the H-CF state as 
a function of SVi in the same system: solid line-lowest in 
doublet; dashed line-next lowest. 

a function of 8V\ (8V3 — 0), for two values of w = 2b. 
The same is also shown for another model of non-zero 
thickness effect, which describes a quantum well in the 
X3 direction, as used in some experimental samples; for 
this the thickness is w = 2.75. In addition, we have 
plotted the squared-overlap of each of the two states in 
the doublet with the single trial ground state for the H- 
CF phase, for w = 2. For the latter, we have used the 
numerically-obtained ground state for a very large posi- 
tive value of 8V\. In that regime, which is similar to the 
LLL Coulomb pseudopotential values, the ground state 
is believed to be incompressible and to lie in the H-CF 
phase. The H-CF state has exactly zero overlap with one 
member of the doublet, while the overlap of the other 
with the H-CF state is not zero. The overlap of the H-CF 
state with the lowest K = state drops abruptly to zero 
for 8V1 less than about 0, which is due to a level cross- 
ing. Beyond that point, the other K = state is lower, 
and has zero overlap with the H-CF state, but the over- 
lap of the other with the H-CF state remains non-zero 
and is shown as the dashed line. This behavior, both the 
vanishing overlap and the level crossing, indicates that 
the two members of the doublet have different rotational 
symmetry in the hexagonal geometry, and only one has 
the same symmetry as the H-CF ground state. The non- 
zero overlap with the H-CF state is much smaller than 
that with the two-dimensional RR trial subspace, but in- 
creases with increasing r5Vi, while the overlaps with RR 
slowly decrease (this behavior is consistent with that ob- 
served on the sphere at smaller sizes |8j ; note that on the 
sphere the RR and H-CF states occur at different N^, 



FIG. 3: The low-lying energy spectrum as a function of SVi 
on the torus, for 15 particles. The K = and K / levels 
are shown as distinct symbols. 

which is not the case on the torus). These results suggest 
that there must be a phase transition in the thermody- 
namic limit between these two incompressible phases at 
some value of 8V\, and this occurs for a range of values 
of thickness w. 

In Fig. |3l we show the low-lying energy spectrum for 
N = 15 particles as the 8V\ is varied away zero, for 
w = 2, and 8V3 — 0. The level-crossing of the lowest 
two K = levels, mentioned in the previous paragraph, 
can be seen at 8V\ close to 0. Over the approximate 
range —0.01 < 8V1 < 0, the lowest two K = states 
stay remarkably close in energy, and separated from all 
other states in the system. For 8V1 negative, a tran- 
sition occurs at around —0.01 (a similar transition was 
seen in Refs. Beyond that point, the spectrum 

as a function of aspect ratio of the torus (not shown) 
shows signs that the system is in a stripe phase. Fig. 
13 reveals that this transition is the strongest feature in 
the spectrum, which dominates the behavior of moder- 
ately low-lying levels even rather far away in 8V\. The 
behavior of the spectrum, with many levels converging 
to zero at the same point, suggests that this transition 
may be second-order in the thermodynamic limit. In 
addition, it indicates that the energy scales in the RR 
phase are smaller than in the H-CF phase which occurs 
at large positive 8V\. The transition between the latter 
two phases (where the higher K = state moves up into 
the continuum) occurs not far from the level crossing of 
ground states. In contrast to the rapid rise in energy of 
one of the two ground states of the RR region on enter- 
ing the H-CF region, the higher energy levels show only 
gradual changes in this range of 8V1 . This might suggest 
that this transition is first-order in the thermodynamic 
limit, but see also the spectra for 18 particles below. 

In Fig. 0] we show a 3-dimensional plot of the total 
squared-overlap with the RR trial subspace for TV = 15 
particles on the torus, with w = 2, as a function of both 
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FIG. 4: A 3-dimensional plot of the total squared-overlap 
with the RR trial subspace for N = 15 particles on the torus, 
with w — 2, as a function of both 5Vi and SVz- 



FIG. 6: A contour plot of the level splitting between the two 
lowest K = states for 15 particles as a function of SVi and 
SV 3 . 
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FIG. 5: The same as Fig. 2] but shown as a contour plot. 

5Vi and 6V3. This shows that the high, broad maximum 
as a function of 5V± persists over a range of 5Vs , though 
the location shifts. The same information is shown in Fig. 
as a contour plot. We note that the squared-overlap 
falls rapidly towards zero as the stripe region is entered. 
In view of the large dimensions of the Hilbert spaces 
of K = states, the largest squared-overlaps (> 1.40) 
should be viewed as significantly large. 

In Fig. we show the splitting between the two lowest 
K = states for 15 particles as a function of 8V\ and SV3, 
also for w — 2. For 8V3 zero or negative, the splitting is 
smallest close to where the overlap with the RR states 
is largest, however at larger SV3 the region of smallest 
splitting is seen to bifurcate, while the overlaps there are 
unaffected. 

Finally, we show in Fig. [7| the low- lying spectrum for 
N = 18 particles as a function of SVi, for 8Vz = 0. Over- 
all features are similar to those for 15 particles in Fig. 
13 however, here the splitting of the lowest two K = 
states in the possible RR region is not as small as for 
N = 15, and there are now two crossings between these 



FIG. 7: Low-lying energy spectrum, shown as the difference 
from the ground state energy, for 18 particles as a function 
of 8V1. Again, states at K = and K / are shown with 
distinct symbols. 

two states between the strip region and the H-CF region 
(the spectrum in Fig. ^ is taken near one of these cross- 
ings). In addition, the gap to the non-zero K levels de- 
creases somewhat in the transition region to the H-CF, 
and is even comparable to the splitting of the doublet 
in that region (however the gap to other K = states 
remains larger). Nonetheless, we expect that overlaps 
with the RR states remain large, and those with the H- 
CF state, small, in the central region with a strong gap 
above the ground state. We note that the appearance 
of two level crossings, with a similar not-so-small split- 
ting between them also seems to occur in the N = 15 
case if we examine larger values of SV3 (the "bifurcation" 
noted above), while overlaps with the RR subspace re- 
main large. Consequently, this aspect of the N = 18 data 
seems consistent with behavior at N = 15. 

Uniting all aspects of our results, we conclude that 
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over a range of parameter values that includes the n = 
1 Coulomb interaction and experimentally relevant non- 
zero thickness effects, there is large overlap of the lowest 
two K = states with the trial subspace, and relatively 
small splittings of those two states, suggesting that they 
become degenerate in the thermodynamic limit. This 
region is bordered on one side by the H-CF phase, which 
has poor overlap (and for some parameters zero overlap, 
due to different symmetry) with the lowest member of 
the doublet in this region. It is bordered on the other 
side by a stripe phase, similar to others in the higher 
LLs. The overall pattern of behavior is quite similar to 
that observed for the MR phase at 5/2 and is also 
consistent with our previous results on the present system 
on the sphere 

To conclude, we find the evidence that the observed 
12/5 state is in the RR phase compelling. The non- 
observation of the 13/5 state so far suggests that LL 
mixing plays a role at the latter filling factor. While a 
systematic study of energy gaps for charged excitations 
will have to await larger system sizes, we find indications 
that the scales in this state will be smaller than in the 
H-CF phase at the same density of particles, which is 
broadly consistent with experimental findings 
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